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1.  Introduction 


Consider  an  experimental  situation  where  it  is  desired  to  compare 
p  >,  2  test  treatments  to  a  control  treatment.  Let  the  p  +  1  treatments 
be  indexed  0,1,..., p  with  0  denoting  the  control  treatment  and  l,2,...,p 
denoting  the  test  treatments.  It  is  desired  to  compare  simultaneously  the 
p  test  treatments  to  the  control.  For  improving  the  precision  of  the  com¬ 
parisons  the  experimental  units  are  to  be  blocked  in  b  blocks  each  of  size 


2  _<  k  ^  p.  We  are  then  in  an  incomplete  block  design  setting. 

Let  Y...  denote  the  observation  on  treatment  i  (0  <  i  <  p) 

i  j  h  —  — 

in  block  j(l  <  j  <  b)  in  plot  h(  1  <.  h  <.  k).  We  assume  the  usual  additive 
linear  model  without  interactions,  namely 


(l.'O 


=  u  +  a.  +  a.  +  e..., 
M  i  pj  i  jh  * 


the  are  assumed  to  be  uncorrelated  random  variables  with  mean  0  and 
2 

common  variance  a  .  The  p  control -treatment  contrasts  -  a.  are  to  be 
estimated  by  their  BLUEs  5q  -  (1  <_  i  <n  p).  It  is  desired  to  choose 

an  experimental  design  (an  allocation  of  treatments  to  blocks)  which  will 
yield  the  best,  in  some  sense,  set  of  estimates  among  all  possible  designs. 

For  given  values  of  b,  k,  and  p  let  C(b,k,p)  denote  the  class  of  all 
possible  incomplete  block  designs  with  b  blocks,  each  of  size  k(p  >  k  >  2), 
p  test  treatments  indexed  l,...,p,  and  a  control  treatment  indexed  0. 


For  a  design  d€C(b,k,p)  let  r. .(d)  denote  the  number  of  replications  of 

J  b 

treatment  i(o  <  i  <_  p)  in  block  j(l  <  j  <_  b).  Also  let  r.(d)  =  [  r.-(d) 

1  i=l  1J 

b  J 

and  x.„(d)  *  f  r^(d)r,.(d)  (0  <  i  t  i  <  p).  Notice  r^(d)  represents 

the  number  of  replications  of  treatments  in  the  entire  design  d  and 


x ^ ^ ( d )  represents  the  number  of  times  treatments  i  and  i  are  paired  to 
gether  in  a  block  summed  over  all  blocks. 
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For  d€C(b,k,p),  let  M(d)  denote  the  information  matrix  correspond¬ 
ing  to  estimating  all  og  -  o^,  1  <  i  <  p,  (as  in  Bechhofer  and  Tamhane 
(1981)).  M(d)  is  a  nonnegative  definite  p*p  matrix  and  is  nonsingular 
if  and  only  if  all  the  ag  -  a^  are  estimable,  in  which  case  it  is  proportional 
to  the  inverse  of  the  covariance  matrix  of  ag  -  c^,  1  <_  i  <_  p. 

We  now  make  our  goal  of  finding  a  design  d€C(b,k,p)  which  gives  us 

the  best,  in  some  sense,  set  of  BLUEs  ag  -  a^(l  <_  i  <  p)  more  explicit. 

Following  the  work  of  Kiefer  (see  for  example  Kiefer  (1958,  1959,  1971, 

and  1974)) we  seek  a  d€C(b,k,p)  which  minimizes  <j>(M(d))  for  some  function 

<P  over  C(b,k,p).  Such  a  design  will  be  called  $ -optimal.  Restricting  to 

non-singular  designs,  some  common  examples  of  <t>  are  $g(M(d))  *  detM~^(d) 

(so  called  D-optimality) ,  <|>,(M(d))«  tr  M"^(d)  (so  called  A-optimal ity) , 

and  4>oo(M(d))  *  maximum  eigenvalue  of  M(d)  (so  called  E-optimality).  In 

the  present  context  of  control -treatment  comparisons,  A-optimality  has  an 

P  .  A 

appealing  statistical  interpretation,  viz.  it  minimizes  ^  varfag-a^) 

over  all  designs.  We  are,  however,  yet  to  realize  natural  statistical 
interpretations  for  the  other  criteria. 

Traditionally,  Kiefer  and  other  researchers  were  interested  in  an 
orthonormal  basis  of  treatment  contrasts.  In  other  words,  the  aim  was  to 
determine  good  designs  for  estimating  Pa  where  a  is  the  vector  of  all  the 
p  +  1  treatment  effects  and  P  is  a  p  *  p  +  1  matrix  of  zero  row  sums  and 
orthonormal  rows.  Nothing  much  seems  to  be  known  for  the  situation  when 
the  contrasts  are  not  mutually  orthogonal.  In  this  paper  we  look  at  one 
such  situation  -  that  of  control -treatment  comparison. 

Let  y.(d),  1  <  1  <_  p,  be  the  positive  eigenvalues  of  the  well  known 
"C-matrix"  of  normal  equations  for  a,  for  a  design  d  in  p  +  1  treatments 
in  b  blocks  of  k  plots  each.  Let  Pa  be  any  vector  of  p  independent 
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treatment  contrasts,  and  V(Pa(d))  be  the  covariance  matrix  of  the  BLUE's 
of  Pa.  Then  it  can  be  shown  that 

det  V(Pa(d))  =  (det(PP,))(v1(d)...vp(d))‘1. 

This  can  be  established  by  starting  from  a  spectral  decomposition  of  the 
C-matrix  for  the  design  d,  or  by  proving  a  result  like  equation  (A. 2)  of 
Bechhofer  and  Tamhane  (1981).  Since  a  B.I.B.  design,  if  it  exists,  is 
D-optimal  in  the  traditional  sense  of  estimating  orthonormal  contrasts, 
we  have  the  following  theorem. 

Theorem  1 .1  A  B.I.B.  design,  if  it  exists,  is  D-optimal  for  estimating 
any  set  of  p  independent  treatment  contrasts. 

It  has  come  to  our  notice  that  this  result  has  been  known  for  some  time 
(Hedayat  (1974)).  Observe  that  the  D-optimal ity  criterion  ignores  the 
particular  interests  of  the  experimenter  expressed  through  the  matrix  P. 

From  the  work  of  Kiefer  and  others  on  optimal  incomplete  block  de¬ 
signs  for  estimating  an  orthonormal  basis  of  treatment  contrasts,  it  is 
known  that  the  B.I.B.  design  is  optimal,  not  only  according  to  the  D- 
criterion,  but  under  a  very  large  class  of  optimality  criteria  as  well 
(see  Kiefer  (1958,  1959,  1971,  1974  and  1975}).  Such  results  might  lead 
us  to  expect  that  in  our  setting  an  optimal  design  d  in  C(b,k,p)  would 
be  symmetric  (in  some  sense)  and  binary  in  the  test  treatments  l,...,p 
(but  not  in  the  control).  Since  the  control  plays  a  special  role  in  our 
setting  we  might  also  expect  that  the  number  of  replications  of  the  con¬ 
trol  (more  specifically  the  r^d))  will  be  an  important  factor  in  deter¬ 
mining  what  design  d  is  optimal.  These  expectations  are  indeed  found  to 
be  the  case  as  will  be  seen  in  the  results  of  section  2. 
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The  proper  sense  of  symmetry  in  a  design  d€C(b,k,p)  turns  out  to  be 
that  all  x^(d)  are  equal  for  1  1  i  t  £  <_  p  and  all  ^(d)  for  1  <  i  <  p 
are  equal  {but  not  necessarily  to  the  x^d)  for  1  <_  i  <_  p).  Such  designs 
are  called  balanced  treatment  incomplete  block  designs  (abbreviated  BTIBs) 
and  were  first  introduced  in  Bechhofer  and  Tamhane  (1981)  in  connection 
with  making  joint  confidence  statements  about  the  contrasts  aQ  -  c^. , 

1  <  i  <  p.  The  interested  reader  is  referred  to  this  paper  for  more  infor¬ 
mation  on  BTIB  designs.  We  remark  that  if  a  design  d£C(b,k,p)  is  a  BTIB  design 
its  information  matrix  M( d)  is  completely  symmetric  (i.e.  all  off  diagonal 
elements  equal  and  all  diagonal  elements  equal).  Bechhofer  and  Tamhane 
(1981)  also  have  a  review  on  available  literature  for  designs  for  control- 
treatment  comparisons. 

Section  2  of  this  paper  contains  results  about  what  designs  are  «j»- 
optimal  for  a  fairly  broad  class  of  functions  As  an  important  applica¬ 
tion  we  discuss  A-optimal  designs. 

The  class  of  functions  considered  in  section  2  does  not  include 
E-optimality.  This  is  treated  in  section  3,  which  also  includes  a  result 
showing  that  an  A-optimal  design  is  optimal  according  to  another  statistic¬ 
ally  interesting  criterion.  Section  4  contains  some  concluding  remarks. 

2.  A-Qptimal  Designs 

We  begin  this  section  with  a  series  of  lemmas  culminating  in  a  general 
theorem  from  which  A-optimal  designs  may  be  obtained  as  a  special  case. 

Suppose  d€C(b,k,p)  is  arbitrary.  Let  %  be  the  set  of  all  p!  permu¬ 
tations  of  the  test  treatments  l,...,p.  Let  ad,o€$,  be  the  design  re¬ 
sulting  from  d  by  the  permutation  a  of  the  treatments  in  d.  We  define 

(2.1)  M(d)  =  l  M(cd)/p!  *  l  tt *M(d)u/p! 

a€$  tt€H 

where  n  is  the  set  of  all  p  x  p  permutation  matrices. 
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Lemma  2. 1 .  If  d€C(b,k,p)  then  fi(d)  has  eigenvalues  w-j(d),  u2(d)  *  ... 


=  Wp(d)  with 


^  XQi(d)/k)/p  *  (rQ(d)-  r*.(d)/k}/p 

u2(d)  a  {  \  **.(d)-  f  l  r?.(d)/k-(rQ(d)  -  J  roi  (d)/k)/p>  /  (p-1 ). 
c  1*1  1*1  j-1  u  j*l  UJ 

In  addition  if  d  is  binary  in  test  treatments 

b  ? 

u2(d)  =  {b(k-l)-((k-l)/k)r0(d)-(r0(d)-  \  rgj (d)/k)/p>/(p-l ) 

J  1 

£f.  From  the  appendix  of  Bechhofer  and  Tamhane  (1981),  the  entries  of 


M(d)  are 


r,  (d)  -  l  r\  - ( d)/k  (i,*ij 
’l  j«1  V  1  £ 


m  .  = 

V’2 


■xi  ,  (d)/k 
1 1 1 2 


<VV 


and  the  sum  of  the  entries  in  the  i-th  row  (or  i-th  column)  is  xQl-(d)/k. 
Thus  is  it  straightforward  to  check  that 

(2.2)  R(d)  *  ({  l  r.(d)  -  £  £  rf.(d)/k+  £  X.  ,  /k(p-l)}/p)I 

i=l  1  i=l  j=l  Uiyigip  V2  p 

-  (  I  X.  •  /kp(p-l))J 

1<11<12<p  ^ 1 2  P,P 

where  I  is  the  p  «  p  identity  matrix  and  J  is  the  p  *  p  matrix  all  of 
P  P  »P 

whose  entries  are  +1.  The  first  part  of  the  lemma  now  follows  from  the 
well  known  fact  that  alp  +  bJp  p  has  eigenvalues  a  with  multiplicity  p-  1 
and  a  +  bp  with  multiplicity  1.  The  second  part  involves  essentially 
straightforward  computations  only. 
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Lemma  2.2.  Suppose  ^  is  a  convex  real-valued  possibly  infinite  function 
on  the  set  of  all  p  x  p  non-negative  definite  matrices  and  $  is  invariant 
under  permutations,  i.e.  if  it  is  a  permutation  matrix,  <j> ( tt ' Mir )  * 

Then  for  d€C(b,k,p),  ,j>(ffl{d))  <  $(M(d)). 

pf.  $(M{d))  =  l  <j>(M(cd))/p!  since  $  is  permutation  invariant.  Thus 

by  convexity  $(M(d))  >  $(  T  M(od)/p!)  a  <p(M(d)). 

o  6* 

Lemma  2.3.  Suppose  <j>  is  some  real -valued  possibly  infinite  function  on 
the  set  of  all  non-negative  definite  p  «  p  matrices  with  the  property  that 
if  M  and  N  are  non-negative  definite  p  x  p  matrices  with  eigenvalues 
»•]  1^2  -  -  yp  an<*  -  v2  —  ***  —  vp  respectively  which  satisfy 

Uj  >_  Vi  for  i  *  1 ,. . .  ,p  then  <j>(M}  £  $(N). 

Let  deC(b,k,p)  be  a  design  which  is  not  binary  in  test  treatments. 
Then  there  exists  d*€Clb,k,p)  which  is  binary  in  test  treatments  with 
rQ(d*)  *  rQ(d)  and  which  satisties  <j>(R(d*))  <_4>(R(d)). 
pf.  In  each  block  of  M(d)  replace  any  duplicates  of  test  treatments  by 
test  treatments  not  in  the  block  so  that  each  block  is  binary  in  test 
treatments  (this  is  possible  since  k  £  p) .  Call  the  resulting  design  d*. 

Notice  d*  is  binary  in  test  treatments,  has  rgj(d*)  =  rgj(d)  ^or 

P  P 

1  <  j  <  b,  and  has  V  r.(d*)  =  T  r.(d).  As  a  result  it  is  easy  to  see 
"  "  i=l  1  i=l  1 

p  b  ?  p  b  ? 

l  l  ij(d*)  1  l  l  r*(d). 

i*l  j=l  1J  i=l  j=l  1J 


From  lemma  2.1  it  then  follows  that  the  eigenvalues  of  R(d)  and  R(d*) 
satisfy  ^(d)  =  ^(d*)  and  y2(d)  =  ...  =  up( d)  <  u2(d*)  »  ...  s  wp(d*). 
Hence  by  the  property  of  $  given  in  the  statement  of  the  lemma, 

$(fi(d*))  <  $(fi(d)). 
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Lemma  2.4.  Among  all  non-negative  integers  rQi »r02** • • ,rob  satisfying 

b  b  2 

l  rn.  =  r,  where  r  is  a  fixed  constant,  the  value  of  T  rf.  is  minimized 
j*l  UJ  j=l  00 

by  choosing  r  -  b[r/b]  of  the  rgj  to  have  value  [r/b]  +  1  and  the  remain¬ 
ing  b(l+[r/b])  -  r  of  the  rg.  to  have  value  [r/b].  Here  [•]  denotes  the 
greatest  integer  function. 
pf.  This  is  leinna  2.3  of  Cheng  and  Wu  (1980). 

Lemma  2.5.  Suppose  $  is  as  in  lemma  2.3.  Suppose  d€C(b,k,p)  is  binary 
in  test  treatments  and  has  rg(d)  >  bk/2.  Then  there  exists  d*€C(b,k,p) 
which  is  binary  in  test  treatments,  has  rg(d*)  <  bk/2,  and  satisfies 
♦(B(d*))  <  ♦(M(d)). 

pf.  Take  d*  to  be  the  design  where  in  each  block  of  d  we  replace  all 
test  treatments  by  the  control  and  all  of  the  original  replications  of  the 
control  by  differing  test  treatments  not  originally  in  the  block.  Notice 
b  b 

(2.3)  rQ(d*)  =  ^  rQj(d*)  =  (k-rQ.(d))  =  bk  -  rQ(d)  <  bk/2  <  rQ(d) 

b  o  b  ^ 

(2.4)  r  (d*)  -  l  r‘.(d*)/k  =  rQ(d)  -  l  r^(d)/k. 

j=l  j=l  J 

From  (2.4)  it  follows  that  if  uj(d),  u2(d)  =  •••  *  up(d)  and 
u2(d*)  *  ...  =  up(d*)  are  the  eigenvalues  of  fl(d)  and  ft(d*),  respectively, 
as  given  in  lemma  2.1,  then  u-|(d)  *  w-j(d*).  Also  from  (2.3),  (2.4),  and 
lemma  2.1,  u2(d)  <  u2(d*).  By  the  property  of  $  given  in  the  lemma  it 
follows  that  $(M(d*))  <_$(M(d)). 

Theorem  2.1.  Suppose  41  is  a  real -valued  possibly  infinite  function  on 
the  set  of  all  p  x  p  non-negative  definite  matrices  satisfying 
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♦(M)  »  l  f(vO 
i=l  1 

where  <.  1  •••  i  up  are  the  eigenvalues  of  M,  f  is  a  real  valued 

possibly  infinite  function  on  the  set  of  all  non-negative  numbers  which 
Is  continuous  on  the  set  of  all  positive  numbers,  has  f'  <  0  and  f"  >  0 
(here  primes  denote  differentiation).  Suppose  there  is  a  5€C(b,k,p)  such 
that  M( 6 )  is  completely  symmetric  and 

(i)  6  is  binary  in  test  treatments 

(ii)  rQ(6)  is  the  value  of  the  integer  r,  0  £  r  £  [bk/2],  which 
minimizes 

(2.5)  g(r;b,k,p)  =  f((r-h(r;b)/k)/p) 

+  (p-1  )f( (b(k-l  )-((k-l )/k)r-(r-h(r;b)/k)/p)/(p-l ) ) 

where 

(2.6)  h(r;b)  *  (b(l+[r/b])-r) [r/b]2  +  (r-b[r/b]) ( [r/b]+l )2 

(iii)  the  ^(6)  have  value  either  [rQ(6)/b]  or  [rQ(6)/b]  +  1. 

Then « is  ^-optimal  over  C(b,k,p). 


pf.  First  we  notice  $  and  f  have  the  following  properties 

(a)  <j>  is  convex  and  orthogonal  invariant  (i.e.  if  ir  is  an  orthogonal 

matrix  then  <<> ( tt ' M-rr )  =  $(M)) 


j,  1 ,!, 


f(v.j) 


Ui  > 


for  all 


(c)  if  p'j  ~ 


V  V1  -  v2  = 


=  v  ,  y-j  >_  v-j ,  and 


J  ^  v1  then  ^  f(ui)  £  ^  f(v.), 


Property  (a)  follows  from  the  fact  that  f"  >  0.  Property  (b)  follows  from 
the  fact  that  f'  <  0.  Property  (c)  follows  from  the  fact  that  <j>  regarded 
as  a  function  of  (^.....u  )'  is  Schur  convex. 
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Now  suppose  6  is  as  given  in  the  theorem.  Let  d€C(b,k,p)  be  any 
design  which  is  binary  in  the  test  treatments.  By  (a)  and  lemma  2.2, 

<|>(M(d))  £$(M(d)).  If  rQ(d)  >  bk/2  by  lemma  2.5  there  exists  d*€C(b,k,p) 
with  r*Q(d*)  <  bk/2  and  <j>(M{d*))<  $(M(d)).  Replace  d  by  d*.  If  ^(d)  <_  bk/2 
let  d*  *  d.  Notice  by  lemma  2.1  M(d*)  has  eigenvalue? 

b  p 

u-j (d*)  =  (rQ(d*)  -  ^  rgj(d*)/k)/p 

i^d*)  =  ...  =  up(d*) 

b  ? 

■  {b(k-l)-((k-l)/k)r0(d*)-(r0(d*)-  roj(d*)/k)/p}/(p-U. 

Using  the  facts  k  ^  2,  <"g(d*)  <_  bk/2,  and  some  calculus,  one  can  prove 

that. 


y2(d*)  >  bk/4p 


2 

=  max(rQ-  J  rQ./k)/p 
j=l  J 

1  u-j  (d*) 


where  the  maximum  is  over  all  real  numbers  rg*rg] »• • • *rgj  such  that 

r0  i  0*  rQj  1  1  1  J  i  b»  anb  1  roj  =  ro  -  bk/2*  Tbus  tbe  eigenvalues 

j  i 

u-,  (d*)  ,n2(d*) ,. . .  ,up(d*)  of  fi(d*)  satisfy  Ul(d*)  <  u2(d*)  *  ...  =  up(d*). 

Next  notice 
(2.7)  <j>(M(d*) ) 

b  9 

=  f(Ul(d*)>  +  (p-l)f(u2(d*))  =  f({ rQ(d*)  -  l  roj(d*)/k)/P) 

J  ^ 


♦  (p-l)f({b(k-l)-((k-l)/k)rQ(d*)-(rQ(d*)  -  l  r*  (d*)/k)/p>/(p-l )) . 

j=l  J 


For  a  fixed  value  of  rQ( d*)  <_  bk/2  we  have  i^(d*)  +  (p-ljugfd*) 

=  (k-1  )(b-rQ(d*)/k)  =  constant  and  the  largest  possible  value  of  w-j ( d*) 

b  2 

*  (rQ(d*)-  l  rgj(d*)/k)/p  occurs  when  b(l+[r0(d*)/b])-r0(d*)  of  the 

rQj(d*)  are  [rQ(d*)/b3  and  rQ(d*)  -  b[rQ(d*)/b]  of  the  rQj(d*)  are 
[rg(d*)/b]  +  1  by  lemma  2.4.  This  choice  of  the  rgj(d*)  maximizes  uj(d*) 
for  fixed  rQ(d*)  and  hence  by  property  (c)  of  <j>  minimizes  <t>(fi(d*)).  If 
we  then  select  a  value  of  rg(d*)  <.  bk/2  (with  the  optimal  choice  of  the 
rQj(d*))  which  minimizes  the  R.H.S.  of  (2.7)  we  see  that  this  is  precise¬ 
ly  the  value  of  rQ(5)  and  the  rg j ( 5 )  stated  in  the  theorem.  We  thus  con¬ 
clude  6  is  a  design  minimizing  #(fl(d*))  among  all  d*  which  are  binary  in 
test  treatments  and  have  rg(d*)  <  bk/2.  Since  M(<$)  is  completely  symmet¬ 
ric,  M( a )  =  M( 6 )  and  we  see  using  lemma  2.2  that  <j>(M(6 ) )  _<  $(M(d))  <_  4>(M(d)) 
(d  is  the  design,  binary  in  test  treatments,  we  chose  arbitrarily)  we  con¬ 
clude  6  is  ^-optimal  among  all  d  which  are  binary  in  test  treatments. 
Property  (b)  of  and  lemma  2.3  then  give  us  that  6  is  ^-optimal  among  all 
designs. 

Theorem  2.1  is  useful  for  finding  optimal  designs  for  many  <j>  such 

P  P 

as  <j>(M)  =  l  -  Inp.  (D -optimal ity)  and  ij>(M)  =  £  1/y.  (A-optimal ity) . 
i=l  i=l 

It  is  not  directly  applicable  to  the  problem  of  finding  E-optimal  designs, 
i.e.  the  design  $€C(b,k,p)  which  minimizes  the  maximum  eigenvalue  of 
M"^(5)  (or  maximizes  the  minimum  eigenvalue  of  M(<$)). 

As  mentioned  in  the  introduction  A-optimal  designs  are  statistically 
very  meaningful.  So  we  examine  such  designs  in  some  detail.  A  design 
d€C(b,k,p)  is  A-optimal  if  it  minimizes  tr  M"^(d)  over  C(b,k,p).  In  the 
notation  of  theorem  2.1  this  means  <*>(M(d))  *  tr  M"^(d)  and  f(y)  *  1/u. 
Equations  (2.5)  and  (2.6)  then  become 
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(2.7)  g(r;b,k,p)  *  p/{r-(h(r;b)/k)} 

+  (p-l)2/(b(k-l)-r(k-l)/k-(r-h(r;b)/k)/p} 

with 

(2.8)  h(r;b)  *  [r/b]2(b+b[r/b]-r)  +  (r-b[r/b])([r/b]+l)2. 

The  following  result  result  is  a  consequence  of  theorem  2.1. 

Theorem  2.2.  Suppose  R  Is  the  value  of  the  integer  r,  0  £  r  <_  [bk/2], 
which  minimizes  g(r;b,k,p)  as  given  in  (2.7).  Also  suppose  6€C(b,k,p) 
is  a  B.T.I.B.  design  such  that 

(i)  5  is  binary  in  test  treatments 

(ii)  rQ(«)  a  R 

(iii)  rgj ( 6 )  =  [R/b]  or  [R/b]  +  1  for  1  <_  j  <  p 
then  5  is  A-optimal  over  C(b,k,p). 

The  integer  r  which  minimizes  g(r;b,k,p)  can  easily  be  found  using 
a  computer.  As  an  example,  r  =  18  minimizes  g(r;24,3,9)  and  the  follow¬ 
ing  B.T.I.B.  design  is  therefore  A-optimal. 

/0  000000000000000001  1  234  5  \ 

1  1  1  1  222233344566772634581 

V3  4584578579696898997687  9/ 

Here  columns  correspond  to  blocks  and  the  numbers  are  the  treatment  labels. 

Having  determined  the  integer  R,  the  next  step  is  to  investigate 

whether  a  B.T.I.B.  design  satisfying  (i)-(iii)  exists  or  not.  Writing 

q  =  [R/b]  and  a  *  R  -  bq,  an  A-optimal  design  looks  like 


where  d^  consists  of  q  plots  in  each  of  b  blocks  and  d^  the  rest  of  the 
k  -  q  plots  in  the  blocks,  d^  consists  entirely  of  the  control,  while  d^ 
is  binary  in  all  p  +  1  treatments  with  the  control  appearing  a  times. 
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The  A-optimal  design  shown  above  gives  an  example  of  d  *  d(2)  since  here 

q  -  0. 

If  a  =  0,  then  d^)  has  to  be  a  B.I.B.  design  in  the  p  test  treatments. 
The  following  table  gives  some  examples  of  A-optimal  designs  having  this 
structure. 

b  k  b  q 

10  3  5  1 

14  4  7  1 

30  2  4  1 

30  4  5  1 

30  4  6  1 

30  5  10  1 

30  6  25  1 

Let  us  denote  by  d^  a  design  in  C(b,k,p)  which  is  a  B.I.B.  in  the 
k  -  q  plots  of  b  blocks  in  the  p  test  treatments,  augmented  by  the  con¬ 
trol  in  each  of  the  remaining  q  plots  of  b  blocks.  Designs  of  this  type 
have  been  mentioned  briefly  by  Cox  (1958,  p.  238);  Pesek  (1974)  has 
look  at  dj.  Neither  of  them  have  considered  these  as  optimal  designs.  The 

interested  reader  may  find  their  efforts  put  in  perspective  in  Bechhofer 
and  Tamhane  (1981).  We  shall  now  show  that  for  many  q,  d^  cannot  be  a 
very  bad  design  -  it  is  at  least  A  -  better  than  a  B.I.B.  design  in  all 
p  +  1  treatments. 

A  B.I.B.  design  is  a  binary  B.T.I.B.  design  with  r^(d)  *  bk/(p+l). 
Moreover,  for  any  B.T.I.B.  design,  tr  M(d)"^  *  g(rg(d) ;b,k,p) .  Hence  we 
look  at  the  sign  of  the  function, 

(2.9)  g-|  (q)  s  g(qb;b,k,p)  -  g(bk/(p+l )  ;b,k,p) 

■  p/{r-(bq2)/k} 

+  (p-l)2/{b(k-1)-r(k-l)/k-(r-(bq2)/k)/p}  -  2p2/{b(k-l)}  , 
with  q  allowed  to  be  positive  integers  only.  If  one  allows  q  to  be  any 
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real  number,  then  It  is  easy  to  see  that  the  polynomial  equation 

g^q)  *  o 

has  no  roots  in  the  interval  [l,(k-l)/2]  of  q.  Moreover  g^(l)  <  0  and 
g1((k-l)/2)  <  0,  but  g-j { k/2)  >  0.  Thus  in  particular  (q)  <  0,  q  *  1,2, 

. . . ,[(k-l )/2].  We  summarize  this  in  the  following  theorem. 

Theorem  2.3.  d^  is  A-better  than  a  B.I.B.  design  in  all  p  +  1  treatments 
for  all  q  *  1 ,2,. . . ,[(k-1 )/2],  whenever  they  exist. 

3.  E-Optima 1  Designs 

We  now  determine  E-optimal  designs. 

Theorem  3.1 .  If  there  exists  6  €C(b,k,p)  such  that 

(i)  every  block  contains  exactly  k/2  replications  of  the  control,  if 
k  is  even,  or  either  [k/2]  or  [k/2]  +  1  replications  of  the  con¬ 
trol  if  k  is  odd 

and 

(ii)  *  ^02^)  *  •••  *  *0p(5) 

then  6  is  E-optimal  over  C(b,k,p). 

pf.  We  first  show  that  xQ1(<5)/k  Is  the  minimum  eigenvalue  of  M( 6 ) .  To 
see  this,  notice  that  the  sum  of  the  entires  in  the  i-th  row  of  M< <5 )  is 
x01(5)/k.  Since  Xq-j (6)  ■  xoi(s)  for  a11  the  row  sums  of  are 

Xgi(6)/k.  It  therefore  follows  that  the  p  *  1  vector  (1,1, . . . ,  1 ) '  is  an 
eigenvector  of  M (5)  with  eigenvalue  Ag^(<5)/k. 

To  verify  that  xQ1(6)/k  is  the  smallest  eigenvalue  of  M(5),  let  e 
be  any  eigenvector  of  M(5)  other  than  '.  Without  loss  of  gen¬ 

erality  we  may  assume  the  largest  coordinate  in  absolute  value  of  e  is 
+1 .  Suppose  +1  is  the  i-th  coordinate  of  e.  Let  x  denote  the  eigenvalue 


of  M(5)  corresponding  to  the  eigenvector  e.  Let  e-  denote  the  j-th  co- 

v 

ordinate  of  e  and  j ( 6)  the  i,j-th  entry  of  M(6).  The  i-th  coordinate 
of  H(5)  e  is 

(«)  ♦  f  e^m^ . ( 6)  >  +  f  m . . ( 5 ) 

j=l  3  11  j=l  1J 

jfi  jfi 

*  I  mi j ( <5)  =  ^0i(5)/k  *  X01(6)/k. 

J  ^ 

The  Inequalities  above  follow  from  the  fact  that  e*  *  +1 ,  |e.|  £  1  for 

•  j 

all  j,  and  1^.(5)  _<  0  for  i  t  j.  Since  xe  ■  M($)e  and  the  i-th  coordinate 
of  xe  is  xe^  *  x,  we  have  x  _>  xQ1  ( d)/k. 

Since  e  was  an  arbitrary  eigenvector,  hence  x  was  an  arbitrary  eigen¬ 
value,  we  conclude  Xq^(6)/Ic  is  the  minimum  eigenvalue  of  M(6). 

Notice 


eimii 


A01  ( 5 ) / k  ®  {  ^01(5)/k>/ps  (r0(«)-i  r0j(<s)/k>/p  . 

jBl 


By  a  proof  similar  to  that  used  in  lemma  2.4  one  can  show  that  among  all 
integers  rgi  »ro2”  *  •  ,rob  sucl1  that  0  ~  r0j  -  the  value  of 
b  2  b 

r0  ”  r0j/k(r0*  ^1  r0j ^  max'i,n*zecl  by  choosing  all  rQj  *  k/2  if  k 

is  even  or  rQj  »  [k/2]  or  [k/2]  +  1,  j  *  l,...,b,  if  k  is  odd.  Since  these 
are  precisely  the  values  of  the  ^(5)  we  conclude 


(3.1)  min  eigenvalue  of  H( 5)  *  (rn(«)  -  T  rn.U)/k)/p 

0  jsj  OJ 

b  o 

l(r0(d)-  l  r£.(d)/k)/p 
j»l  J 


for  all  d  €  C(b,k,p) .  For  any  d€C(b,k,p) 


min  eigenvalue  of  M(d) 


■  min  u'M(d)u 

llulM 

<  (1.1 . l)H(d)(l.l,...,l)7p 

■  l  ?  ^^(dj/p  *  l  xQi(d)/kp 

1«1  j«1  i»l  01 

<  min  eigenvalue  of  M(<s) 

where  we  have  used  the  facts  that  V  m..(d)  »  xn4(d)/k, 

j=l  01 

XQi(d)/k  *  rQ(d)  -  I  ^  rgj (d) ,  and  (3.1). 

We  conclude  that  6  has  the  largest  minimum  eigenvalue  among  all 
d€C(b,k,p)  and  hence  It  follows  that  6  Is  E-optimal  over  C(b,k,p). 

Theorem  3.1  can  be  used  to  find  E-optlma 1  designs  for  comparing  test 
treatments  to  a  control.  When  k  is  even  each  block  In  an  E-optimal  design 
d  must  have  exactly  k/2  replications  of  the  control.  It  is  straightfor¬ 
ward  to  see  that  a  necessary  and  sufficient  condition  for  a  design  d  to 
be  E-optimal  when  k  is  even  is  that  m  *  bk/2p  be  an  integer,  r^(d)  *  m 
for  all  1  <.  1  <.  P*  and  each  block  contain  exactly  k/2  replications  of  the 
control.  Thus  each  of  the  following  designs  are  E-optimal  when  b«k*p*4 

/O  0  0  0  \  /o  0  0  o\  /o  0  0  0  \ 

0  0  0  0  \  0  0  0  0  \  I  0  0  0  0  \ 

1  1  2  3  4  I  l  1  2  3  4  I  I  1  1  3  3  I 

V  2  3  4  /  \  2  3  4  1  /  \  2  2  4  4  / 

When  k  is  odd  there  Is  more  flexibility  in  finding  designs  which  are 
E-optimal.  In  particular  an  E-optimal  design  d  can  have  either  [k/2]  or 
[k/2]  +  1  replications  of  the  control  in  a  given  block  and  this  flexibil¬ 
ity  allows  one  to  meet  the  condition  \^(d)  *  ...  *  Xgp(d)  TOre  eas^y 
than  when  k  is  even.  For  example,  when  b*  5,  k*  5,  and  p*  6  it  is  easy 


to  check  that  the  condition  xQ1(d)  =  ...  =  x06(d)  cannot  be  met  If  all 
blocks  have  exactly  2  replications  of  the  control  or  if  all  blocks  have 
exactly  three  replications  of  the  control.  However  the  design 

0  0  0  0  0 


has  ^Qi ( d)  s  5  for  all  1  £  1  £  6  and  is  E-optimal. 

One  can  show  that  when  k  is  odd  if  b  *  y  +  w,  where  a  and  e  are 
integers  ^0  such  that  s  *  y[k/2]/p  and  t  a  o>([k/2]+l  )/p  are  both  non¬ 
negative  integers,  then  an  E-optimal  design  d  exists.  In  fact  any  design 
consisting  of  y  blocks  in  which  each  test  treatment  appears  exactly  s 
times  and  the  control  appears  exactly  [k/2]  +  1  times  in  each  of  the  y 
blocks,  and  consisting  of  «  blocks  in  which  each  test  treatment  appears 
exactly  t  times  and  the  control  appears  exactly  [k/2]  times  in  each  of 
the  <a  blocks,  is  E-optimal. 

It  should  also  be  remarked  that  If  an  E-optimal  design  also  happens 
to  be  a  B.T.I.B.  which  Is  binary  in  test  treatments,  then  the  design  may 
have  additional  optimality  properties.  For  example,  the  design 


d  = 


with  b*  3,  k*  3,  pa  3  turns  out  to  be  both  E-optimal  and  A-optlmal. 

Even  though  E-optimality  is  an  effective  way  of  minimizing  a  norm  of 
M(d)’1,  in  the  present  context  it  does  not  seem  to  possess  a  very  natural 
statistical  meaning.  If  one  favours  a  minimax  style  approach  then  possibly 


one  way  Is  to  minimize  the  maximum  variance  of  (oq-o^),  the  maximum  being 
over  1  <  i  <  p,  minimum  over  all  d€C(b,k,p).  In  other  words,  one  considers 
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the  criterion  <*>( M( d) )  *  maximum  diagonal  entry  of  M"^(d).  Since  $  is  con¬ 
vex  and  permutation  invariant  lemma  2.2  says  that  $(R(d))  <_  $(M(d))  for 
any  d€C(b,k,p).  Furthermore,  since  R(d)  is  completely  symmetric,  it  is 
easily  verified  that  $(R(d))  =  tr  R(d)’Vp*  From  these  two  observations 
it  follows  that  if  <5€C(b,k,p)  is  as  in  theorem  2.2  than  it  is  ^-optimal 
as  well  as  A-optimal.  This  lends  additional  significance  to  the  A-optimal 
des i gns . 

A  class  of  criteria  that  are  sometimes  considered  in  optimal  design 
investigations  (see  Kiefer  (1974))  are  the  <t>q  criteria,  0  <  q  <  «  ,  where 

♦„(M(d))  »  I  u:q 

and  u-j  £  ...  <_  Up  are  the  eigenvalues  of  M(d).  D-optimality  and  E-optimality 

are  limiting  cases  of  these  criteria  in  the  sense  that  lim(*  (M(d))/p)^q 

q-0  q 

*  (det  M_1(d))1//p  and  1im(#  (M(d))/p)^q  *  max  eigenvalue  of  M_1(d).  In 

q-H»  ° 

particular,  and  ^  are  sometimes  used  to  denote  the  D-optimality  and 
E-optimality  criteria,  respectively.  Also  notice  that  $-j  is  just  the 
A-optimal ity  criterion. 

An  examination  of  our  above  results  for  D-,  A-,  and  E-optimality  (or 
for  ,  and  ^  )  indicates  that  the  number  of  replications  of  the  con¬ 

trol  in  an  optimal  design  is  smallest  for  D-optimality,  second  smallest 
for  A-optimality,  and  largest  for  E-optimality.  This  suggests  that  the 
number  of  replications  of  the  control  in  a  B.T.I.B.  design  which  is  (frq- 
optimal  is  increasing  as  q  increases.  Since  $  ,  0  <  q  <  »  ,  satisfies 
the  conditions  of  theorem  2.1,  It  is  possible  (although  somewhat  tedious) 
to  verify  that  this  is  indeed  the  case.  From  this  it  follows  that  if 
d€C(b,k,p)  is  a  B.T.I.B.  design  which  is  binary  in  test  treatments  with 
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[blytp+l)]<rg(d)  <  b[k/2],  then  d  Is  ^-optimal  for  some  0  <  q  <  ». 

4.  Concluding  Remarks 

B.T.I.B.  designs  were  introduced  by  Bechhofer  and  Tamhane  (1981). 

The  property  of  being  balanced  in  test  treatments  make  these  designs 
attractive  for  use  in  comparing  test  treatments  to  a  control.  It  has 
been  an  open  question  as  to  what  sort  of  optimality  properties,  if  any, 
these  designs  might  have.  This  paper  answers  that  question  in  part  by 
showing  that  certain  B.T.I.B.  designs  do  indeed  have  some  optimal  pro¬ 
perties.  It  is  hoped  that  this  paper  will  provide  added  incentive  for 
the  study  of  these  designs,  particularly  their  construction. 

One  interesting  feature  of  the  results  contained  in  this  paper  is 
that  exactly  what  design  is  optimal  depends  very  mcuh  on  the  optimality 
criterion  used.  This  is  different  from  the  usual  incomplete  block  de¬ 
sign  setting  where  orthonormal  treatment  contrasts  are  of  interest.  In  that 
setting  optimal  designs  are  often  found  to  be  somewhat  independent  of 
the  criterion  used  and  hence  selection  of  an  appropriate  design  is 
simplified.  The  results  of  this  paper  show  that  the  selection  of  an 
appropriate  design  for  comparing  test  treatments  with  a  control  is  un¬ 
fortunately  sensitive  to  the  optimality  criterion  used.  Unless  one  has 
a  clear  idea  of  the  optimality  criterion  one  wishes  to  use  in  selecting 
a  design,  selection  of  an  appropriate  design  is  not  as  straightforward 
as  in  the  traditional  setting  of  estimating  an  orthonormal  basis  of 
treatment  contrasts.  We  do  recognize,  however,  the  poignant  statistical 
appeal  of  A-optimality. 
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pf.  In  each  block  of  M(d)  replace  any  duplicates  of  test  treatments  by 
test  treatments  not  in  the  block  so  that  each  block  is  binary  in  test 
treatments  (this  is  possible  since  k  ±  p).  Call  the  resulting  design  d*. 
Notice  d*  is  binary  in  test  treatments,  has  ^(d*)  =  rQ j ( d )  for  all 

1  <  j  <  b,  and  has  7  r - ( d*)  =  7  r . ( d ) .  As  a  result  it  is  easy  to  see 
i=l  1  i=l  1 


1=1  J=1  J  1=1  J=1  J 


From  lemma  2.1  it  then  follows  that  the  eigenvalues  of  M(d)  and  fi(d*) 


satisfy  u-j (d)  =  ^(d*)  and  y2(d)  =  ...  =  up(d)  <_  u2(d*)  =  ...  =  up(d*). 
Hence  by  the  property  of  given  in  the  statement  of  the  lemma, 

♦  (A(d*))  <  *(«(<!)). 
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